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Abstract 

(N 

In this paper it is shown that the Brans - Dicke scalar field itself can serve the purpose 
of providing an early deceleration and a late time acceleration of the universe without any 
need of quintessence field if one considers an interaction, i.e, transfer of energy between the 
dark matter and the Brans - Dicke scalar field. 

in 
o 

o : 1 Introduction 



Over the last few years, the speculation that our universe is undergoing an accelerated expansion 
has turned into a conviction. The recent observations regarding the luminosity - redshift relation 
of type la supernovae [1] and also the observations on Cosmic Microwave Background Radiation 
(CMBR) [2] very strongly indicate this acceleration. These observations naturally lead to the 
search for some kind of matter field which would generate sufficient negative pressure to drive 
the present acceleration. Furthermore, observations reveal that this unknown form of matter, 
popularly referred to as the "dark energy" , accounts for almost 70% of the present energy of the 
universe. This is confirmed by the very recent Wilkinson Microwave Anisotropy Probe (WMAP) 
data [3]. A large number of possible candidates for this "dark energy" component has already 
been proposed and their behaviour have been studied extensively. There are excellent reviews on 
this topic [4]. 

It deserves mention that this alleged acceleration should only be a very recent phenomenon 
and the universe must have undergone a deceleration (deceleration parameter q = — J^° 2 > 0) 
in the early phase of matter dominated era. This is crucial for the successful nucleosynthesis as 
well as for the structure formation of the universe. There are observational evidences too that 
beyond a certain value of the redshift z ( z ~ 1.5 ), the universe surely had a decelerated phase 
of expansion [5]. So, the dark energy component should have evolved in such a way that its effect 
on the dynamics of the universe is dominant only during later stages of the matter dominated 
epoch. A recent work by Padmanabhan and Roy Choudhury [6] shows that in view of the error 
bars in the observations, this signature flip in q is essential for the conclusion that the present 
universe is accelerating. 
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So, we are very much in need of some form of a field as the candidate for dark energy, which 
should govern the dynamics of the universe in such a way that the deceleration parameter q was 
positive in the early phases of the matter dominated era and becomes negative during the later 
stages of evolution. One of the favoured choices for the "dark energy" component is a scalar field 
called a quintessence field ( Q-field ) which slowly rolls down its potential such that the potential 
term dominates over the kinetic term and thus generates sufficient negative pressure for driving 
the acceleration. A large number of quintessence potentials have appeared in the literature and 
their behaviour have been studied extensively ( for a comprehensive review, see [7] ). However, 
most of the quintessence potentials do not have a proper physical background explaining their 
genesis. In the absence of a proper theoretical plea for introducing a particular Q-field, non- 
minimally coupled scalar field theories become attractive for carrying out the possible role of the 
driver of the late time acceleration. The reason is simple; the required scalar field is already there 
in the purview of the theory and does not need to be put in by hand. Brans - Dicke theory is 
arguably the most natural choice as the scalar - tensor generalization of general relativity (GR) 
because of its simplicity and a possible reduction to GR in some limit. Obviously Brans - Dicke 
(BD) theory or its modifications have already found some attention as a driver of the present 
cosmic acceleration [8] (see also [9]). It had also been shown that BD theory can potentially 
generate sufficient acceleration in the matter dominated era even without any help from an exotic 
Q - field [10]. But this has problems with the required 'transition' from a decelerated to an 
accelerated phase. Amongst other nonminimally coupled theories, a dilatonic scalar field had also 
been considered as the driver of the present acceleration [11]. 

In most of the models the dark energy and dark matter components are considered to be non- 
interacting and are allowed to evolve independently. However, as the nature of these components 
are not completely known, the interaction between them will indeed provide a more general 
framework to work in. Recently, Zimdahl and Pavon [12] have shown that the interaction between 
dark energy and dark matter can be very useful in solving the coincidence problem ( see also ref 
[13] ). Following this idea, we consider an interaction or 'transfer of energy' between the Brans - 
Dicke scalar field which is a geometrical field and the dark matter. The idea of using a 'transfer' 
of energy between matter and the nonminimally coupled field had been used earlier by Amendola 
[14]. We do it specifically for a modified Brans - Dicke theory. The motivation for introducing 
this modification of Brans - Dicke theory is the following. In the presence of matter and a 
quintessence field, with or without an interaction between them, the evolution of net equation 
of state parameter w plays a crucial role in driving a late surge of accelerated expansion. But 
WMAP survey indicates that the time variation of w may be very severely restricted [15]. If the 
late acceleration is driven by an exchange of energy between matter and a geometrical field <j), 
the question of the variation of w would not arise. 

We write down the Brans - Dicke field equations in the so called Einstein frame. The field 
equations in this version look simpler and G becomes a constant. But one has to sacrifice the 
equivalence principle as the rest mass of a test particle becomes a function of the scalar field 
[16]. So, the geodesic equation is no longer valid and the different physical quantities loose their 
significance. Nevertheless, the equations in this version of the theory enables us to identify the 
energy contributions from different components of matter. However, for final conclusions we go 
back to the original atomic units where we can talk about the features with confidence. We 
choose a particular form of the interaction and show that a constant BD parameter u can not 
give us the required flip from a positive to a negative signature of q in the matter dominated era. 
We attempt to sort out this problem using a modified form of BD theory where u is a function 
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of the scalar field [17]. It has been pointed out by Bartolo and Pietroni [18] that a varying 
uj theory can indeed explain the late time behaviour of the universe. By choosing a particular 
functional form of uj, we show that in the interacting scenario, one can obtain a scale factor 'a' 
in the original version ( i.e, in atomic units ) of the theory so that the deceleration parameter 
q has the desired property of a signature flip without having to invoke any quintessence field in 
the model. We also calculate the statefinder pair {r,s}, recently introduced by Sahni et al [19], 
for this model. The statefinder probes the expansion dynamics of the universe in terms of higher 
derivatives of the scale factor, i.e, a and a. These statefinder parameters along with the SNAP 
data can provide an excellent diagnostic for describing the properties of dark energy component 
in future. 



2 Field Equations and Solutions 

The field equations for a spatially flat Robertson - Walker spacetime in Brans - Dicke theory are 

q a 2 _ p m uj4> 2 d0 

V - + 2 02 6 a(j) ' W 



9- + — --^L-t_i-- (2) 

a a 2 2 2 a ' 



The field equations have been written with the assumption that at the present epoch the 
universe is filled with pressureless dust, i.e, p m = 0. Here p m is the matter density of the universe, 
is the Brans - Dicke scalar field, a is the scale factor of the universe and uj is the BD parameter. 
An overhead dot represents a differentiation with respect to time t. 

The usual matter conservation equation has the form 

p m + 3Hp m = . 

But here we consider an interaction between dark matter and the geometrical scalar field and 
write down the matter conservation equation in the form 

p m + 3Hp m = Q , (3) 

such that the matter field grows or decays at the expense of the BD field. The matter itself is 
not conserved here and the nature of interaction is determined by the functional form of Q. We 
do not use the wave equation for the BD field here because if we treat equations (1), (2) and (3) 
as independent equations, then the wave equation comes out automatically as a consequence of 
the Bianchi identity. It deserves mention that the wave equation will be modified to contain Q 
which will determine the rate of pumping energy from the BD field to matter or vice-versa. This 
interaction term Q is indeed a modification of Brans - Dicke theory. But this interaction does 
not demand any nonminimal coupling between matter and the scalar field and hence does not 
infringe the geodesic equation in anyway. In this interaction, the rest mass of a test particle is 
not modified but rather a "creation" of matter at the expense of the scalar field (or the reverse) 
takes place. In a sense, it has some similarity with the "C - field" of the steady state theory [20]. 

In the Brans - Dicke theory, the effective gravitational constant is given by G = which is 
indeed not a constant. Now, we effect a conformal transformation 
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In the transformed version G becomes a constant. However, this transformation has some limi- 
tations which have been mentioned earlier. But the resulting field equations look more tractable. 
Equations (1) and (2) in the new frame look like 

(4) 

, (5) 

a a 2 4 
and the matter conservation equation takes the form 

Pm + 3 - Pm = Q , (6) 

(X 

where an overbar represents quantities in new frame and ip = Iruf). The scale factor and the 
matter density in the present version are related to those in the original version as 

a 2 = (pa 2 and p m = <p 2 p m . (7) 

Now, we choose the interaction Q of the form 

Q = - a Hp m , (8) 

where a is a positive constant. This negative Q indicates a transfer of energy from the dark- 
matter (DM) component to the geometrical field 0. 

Equation (6) can be easily integrated with the help of equation (8) to yield 

Pm = Po a ( - a - 3) , (9) 

where po is a constant of integration. 

Then, equations (4) and (5) alongwith equation (9) has a solution 

a = At 2 '^ (10) 

where A is a constant given by 

Some arbitrary constants of integration have been put equal to zero while arriving at equation 
(10) for the sake of simplicity. 

Using equations (5) and (10), one can easily arrive at the relation 

.2u) + 3, - 2 4a 1 

If we consider a non- varying u, equation (11) will give rise to a simple power law evolution 
of 0. From equations (7) and (10), the scale factor a in atomic units will also have a power law 
evolution - an ever accelerating or an ever decelerating model contrary to our requirement. This 
is consistent with the exhaustive solutions in Brans - Dicke cosmology obtained by Gurevich et 
al. [21], where the dust solutions are all power law. This indicates that the chioce of constants of 
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integration in equation (10) does not generically change the model. One way out of this problem 
is to consider a generalization of Brans - Dicke theory where the parameter a; is a function of 
the scalar field rather than a constant [17]. An evolving u will be a contributory factor in 
determining the dynamics of the universe. 
We make a choice of uj as, 

2^ + 3 = a , v 

4 ( 3 + a )2( v ^_l)2 • ^ > 

Then, equation (11) can be integrated to yield 

<f> = (1 - fat) 2 , (13) 

0o being a positive constant. 

It deserves mention here that since g 00 and g 00 are both equal to one, the time variable 
transforms as 



This along with equation (13) gives 



__ 1 

00 



dt 2 = (pdt 2 



1 - ^200 (to - t) 



(14) 



which is a monotonically increasing function of t until t = to, beyond which the model really does 
not work. So one can use i itself as the new cosmic time in the original version of the theory 
without any loss of generality. So, for the sake of convenience, from now onwards, we write t in 
place of t. 

We transform the scale factor back to the original units by equation (7), so that we are armed 
with the equivalence principle and can talk about the dynamics quite confidently. 
We have, 

a A t 2 /( 3+ ") 

Also, the Hubble parameter and the deceleration parameter q in the original version comes 
out as, 

J? ?_1 (16) 

3 + a t 1 - o t K ' 

^(l-0ot) 2 -0gt 2 

Q = -1 + 3+a — -5— (17) 

9 + [^(l-0ot)+0ot] 2 • [ 7) 



From equations (15) and (16) it is evident that at t — > ^, both a and H blow up together 
giving a Big Rip. However, this rip has a different characteristic than that engineered by a normal 
phantom field. In the latter, p m goes to zero but pde goes to infinity at the rip. In the present 
case, however, there is no dark energy as such, and the scalar field is a part of geometry and 
hence it is difficult to recognize its contribution to the energy density. In the revised version, 
however, ^^-ib 2 is the contribution towards the stress tensor. It turns out that at t — > -r, this 
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Figure 1: Figure 1(a) and 1(b) shows the plot of q vs. t for different values of a. For figure 1(a) 
we choose a = 0.0000003 whereas for figure 1(b) we set the value as a = 0.0000001. 



contribution remains quite finite. So the big rip is brought into being by the interaction, and not 
by a singularity in the stress tensor. 

The plot of q against t ( figure 1 ) reveals that the deceleration parameter indeed has a sign 
flip in the desired direction and indicates an early deceleration ( q > ) followed by a late time 
acceleration ( q < ) of the universe. Also, the nature of the curve is not crucially sensitive to 
the value of a chosen. 

Equation (5) clearly indicates that a/ a is negative definite. So the signature flip in q = — J^° 2 
in the original Jordan frame must come from the time variation in via equation (15). Local 
astronomical experiments suggest that the present variation of G and hence that of has a very 
stringent upper bound. It is therefore imperative to check whether the present model is consistent 
with that bound. For figure 1, the value of the constant of integration 0o is fixed at 0.3. With 
this value, and the age of the universe taken as ~ 15 Giga years, equation (13) yields 

~ 10~ 10 per year, 

which is consistent with the requirements of the local experiments [22]. The suffix indicates 
the present value. Also in this model, from equations (12) and (15), we get as uo — > oo, — > 1 
and a — > t 2 ^ 3+a K Therefore, for very small value of a (~ 10 -7 ), a is indistinguishable from that 
in GR (a ~ t 2//3 ). This is consistent with the notion that BD theory yields GR in the infinite to 
limit. 



3 Statefinder parameters for the model 

Recently Sahni et al. [19] have introduced a pair of new cosmological parameters {r, s}, termed as 
"statefinder parameters". These parameters can effectively differentiate between different forms 
of dark energy and provide a simple diagnostic regarding whether a particular model fits into the 
basic observational data. These parameters are 



r = J_ an d s = ^_L_ 
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Accordingly, we find the stateflnder parameters for the present model as 

30gt2 _ 3/ 3(! _ ot )2 2(3(1 - cf) tf + 20§t 3 



r = 1 + 



[/3 + o t(l-/^ 



[/3 + o i(l-/3)p 



and 



where (3 = 



3[0^ _ £(1 _ ot )2p + (1 _ /9)0ot] + 2/3 (l _ ot )3 + 



3[/3 + (1 - [-U(3 + (1 - /5)0 o t} + P(l ~ M 2 - 0^ 2 ] 2 



3+a' 



(18) 
(19) 




Figure 2: Plot of r as a function of s for a = .0000001 



If we now plot r(s) for some small value of a ( a << 1 ), we find that the nature of the curve is 
similar to the one expected for scalar field quintessence models with equation of state parameter 
w in the range — 1 < w < [19]. 

4 Discussion 

Thus we see that for a spatially flat universe (k — 0), we can construct a presently accelerating 
model in Brans - Dicke theory or more precisely in a generalized version of it ( uj — oj((f>) ) if 
one considers an interaction between dark-matter and the geometrical scalar field. It deserves 
mention that for this interaction, the Lagrangian should also be modified by the inclusion of an 
interference term between and L m . The salient feature of the model is that no dark energy 
component is required. The nature of the q vs. t curve is not crucially sensitive to small changes 
in the value of a, the parameter which determines the strength of the interaction; only the time 
of 'onset' of acceleration would change by small amounts with a. In revised unit, this interaction 
can be switched off by putting a = 0. However, in Jordan frame it is not possible to switch off 
the interaction with this particular choice. If we put a = 0, the scalar field itself becomes trivial. 
In this frame, the conservation equation ( equation (6) along with equation (8) ) transforms to 

Pm ■ (20) 



p m + 3Hp m = 



aH + 



a — 1 
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As one has both Hp m and |p m in the right hand side, the transfer of energy between matter and 
scalar field takes place both due to the expansion of the universe and the evolution of the scalar 
field. This equation shows that if 

= constant a~ 2a ^ a ^ 

the interaction vanishes. In this case, the transfer of energy due to Hp m and |p m cancel each 
other. 

It is evident from equation (11) that if we consider the interaction between dark matter and 
the geometrical field of the form considered in equation (8), a constant uj will give rise to an ever 
accelerating or ever decelerating model and definitely we are not interested in that. So, the idea 
of varying uj is crucial here as it can very well serve the purpose of providing a signature flip in 
q in this interacting scenario. It deserves mention that the specific choice for the interaction in 
equation (8) and the choice of uj = oj{4>) in equation (12) are taken so as to yield the desired 
result. This is indeed a toy model which simply shows that investigations regarding an interaction 
amongst matter and the nonminimally coupled scalar field is worthwhile. 

From equation (11) it is also evident that (2w + 3) has to be positive definite, i.e, uj has to pick 
up some positive value or at least uj should be greater than — | in order to sustain a consistent 
model. Also, the parameter uj does not have any stringent limit and thus it may be possible to 
adjust the value of uj to some higher value. Equation (12) indicates that if is very close to unity, 
which is consistent with the present value of G, uj can attain a high value at the present epoch. 
This is encouraging as it might be possible to obtain a model which exhibits early deceleration 
and late time acceleration even with a high value of uj, compatible with the limit imposed on it 
by the solar system experiments [22]. 

Also, from equation (13) it is evident that ^ < and | > 0. So, from equation (2), which 
is of particular interest in studying the dynamics of the universe, we see that the last term 
2^| is negative and is the key factor in driving the present acceleration of the universe. This 
term basically provides the effective negative pressure and becomes dominant during the later 
stages of evolution and drives an accelerated expansion. We have also calculated the statefinder 
parameters for the model and show that the {r, s} pair mimics that of a quintessence model. For 
this however, the constant a should be given a very small value (~ 1CT 7 ). In this model we have 
considered a particular form of interaction which indeed is not unique, and some complicated kind 
of interaction may lead to more viable solutions of the various cosmological problems, particularly 
to a model which is not restricted in future. 

Although it is true that General Relativity (GR) is by far the best theory of gravity and the 
natural generalisation of BD theory to GR for large uj limit is shown to be restricted in some 
sense [23], still BD theory always seems to be ready to provide some useful clues to the solution of 
various cosmological problems. The solution to the "graceful exit" problem of inflation in terms of 
an 'extended inflation' scenario [24] was first obtained in BD theory which provided hints towards 
the subsequent resolutions of the problem in GR. Once again here, BD theory in its own right 
could provide a model exhibiting the present cosmic acceleration without introducing any exotic 
dark energy component and since nothing definite is known about the source of this acceleration, 
this type of investigations may lead to some track along which viable solutions may finally be 
arrived at. However, this is a primitive model. This only shows that such investigations can be 
useful. It remains to be seen if the solution is an attractor, and whether the model is consistent 
with the structure formation. 
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